arXiv:1502.06939vl [math.AP] 24 Feb 2015 


Symmetry Breaking and Uniqueness for the 
Incompressible Navier-Stokes Equations 

Radii Dascaliuc * Nicholas Michalowski ' Enrique Thomann* 

Edward C. Waymire^ 


Abstract 

The present article establishes connections between the structure of the deterministic Navier- 
Stokes equations and the structure of (similarity) equations that govern self-similar solutions 
as expected values of certain naturally associated stochastic cascades. A principle result is 
that explosion criteria for the stochastic cascades involved in the probabilistic representations 
of solutions to the respective equations coincide. While the uniqueness problem itself re¬ 
mains unresolved, these connections provide interesting problems and possible methods for 
investigating symmetry breaking and the uniqueness problem for Navier-Stokes equations. In 
particular, new branching Markov chains, including a dilogarithmic branching random walk 
on the multiplicative group (0, oo), naturally arise as a result of this investigation. 

The role of scaling in the question of uniqueness of mild solutions to 3D Navier-Stokes 
equations is the central theme of this investigation. In particular, we describe a framework, 
where the uniqueness for both scale-invariant and general problems is re-cast in terms of 
a non-explosion property of associated stochastic cascades. Thus, if the explosion event of 
the self-similar cascade is probabilistically different from the explosion event for the general, 
non-symmetric cascade in appropriate settings, then we would have a manifestation of sym¬ 
metry breaking in the Navier-Stokes uniqueness problem - the scaling-invariant case being 
qualitatively different. While we are only able to prove partial results related to the associ¬ 
ated explosion problems, the main conclusion of this paper is that the self-similar (scaling- 
invariant) explosion and the general, non-symmetric explosion (in appropriate functional 
settings) are the same, suggesting that scaling symmetry may be directly involved in the even¬ 
tual solution of the outstanding Navier-Stokes well-posedness problem. We note that the idea 
of employing scaling-invariant solutions in the context of well-posedness goes back to Leray 
[25], while the idea to use stochastic cascades to prove existence of mild solutions is due to Le 
Jan and Sznitman [21]. 
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1 Introduction 


1.1 Navier-Stokes equations and their scaling properties. 

The physics of unrestricted three-dimensional incompressible fluid flow is mathematically encoded 
in the corresponding set of Navier-Stokes equations (NSE) governing the time evolution of velocity 
(momentum) u and pressure p in three dimensional Euclidean space. Letting u(x, t) denote the 
velocity of an incompressible fluid at the position x e M 3 and and time t > 0, essentially Newton’s 
law of motion may be cast as 

Qn 

— + u • Vu = z/Au — Vp + g, V • u = 0, u(x, 0 + ) = u 0 (x), x e M 3 , t > 0, (1.1) 

where v > 0 is a positive (viscosity) parameter, V = (O/Axy) i<j< 3 , A = V ■ V is the (vector) 
Laplacian operator, and g is an (external forcing) function with values in M 3 . More generally (1.1) 
may be posed on a domain in R 3 with boundary. However for convenience in this paper we will 
consider the free-space model without boundary or an external force. 

The term du/dt + u ■ Vu represents the acceleration of a fluid parcel within a Lagrangian 
reference frame. In particular, the non-linearity u • Vu is intrinsic to this description of the flow 
and cannot be eliminated. The viscous force z/Au is the result of a linearization of stress-strain 
forces between fluid parcels composing the fluid, and the divergence-free condition V ■ u = 0 
provides conservation of mass; also referred to as incompressibility. The pressure gradient term 
Vp is a fourth unknown in the set of four equations describing the n = 3 coordinates of velocity 
u = (wi, w 2 , w. 3 ) and the scalar pressure p. We refer to [24] and [31] for more background on the 
physical derivation of the Navier-Stokes equations. 

The unique determination of u from the given viscosity parameter v > 0, external forcing g 
(in our case g = 0), and initial data u 0 is an obvious question for both the physics and the mathe¬ 
matics of fluid flow. After more than one-hundred years of research it remains unknown whether 
smooth initial data u 0 leads to the existence of unique smooth (regular) solutions, valid for all 
time. It is believed that mathematical progress on this issue is closely connected to understanding 
the physical phenomenon of turbulence. As a consequence, the resolution of the uniqueness and 
regularity problem for the Navier-Stokes equations ranks among the most important open problems 
of contemporary applied and theoretical mathematics. 

The current state of the regularity issue may be viewed through the prism of natural scaling 
(symmetry) peculiar to of the Navier-Stokes equations as follows 

If u(x, f),p(x, t) is a solution to (1.1), then for any scaling parameter r > 0, 

u r (x, t) = ru(rx, r 2 f),p r (x, t) = r 2 p(rx.,r 2 t ) is also a solution (1.2) 

with initial data ru 0 (rx). 

The quantities of the flow (typically represented by certain norms of u) that preserve this scal¬ 
ing are called critical, the ones that grow as r —* 0 are super-critical, and the ones that decrease 
are sub-critical. For example, since the pioneering work of Leray in the 1930’s (see [25], which 
still remains a benchmark for the regularity problem), it is known that NSE possess global-in-time 
weak solutions that are bounded in L 2 . If we re-scale the L 2 -norm of u according to the scaling 
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above we obtain 


11U-r-11 2,00 ' = sup ||u r (t)||| = sup [ r 2 \u(rx,r 2 t)\ 2 dx=-\\u\\l 00 , 
t£[0,oo) £E[0 ,oo) J T 

R 3 

and so || • || 2 ,oo is a super-critical quantity. Yet, according to Leray’s result, the solution is regular 
as long the L 2 norm of the (vector) gradient remains bounded: 

r 3 ' 

l|Vu r ||| !00 := sup ||Vu r (f)||| 2 = sup / r 2 ^(a xi ^-(rx,r 2 f)) 2 (ix = r||Vu||| jOOJ 

t£[0,oo) tg[0,oo) J ■ ■_| 

R3 

i.e. Leray’s regularity condition is sub-critical. 

More modem regularity criteria still suffer similar scaling defects: The Escauriaza, Seregin and 
Sverak criterion involving boundedness of the L 3 -norm (see [12]), as well as the Koch and Tataru 
condition of smallness of the initial data in the BMO -1 functional space ([20]), are each critical in 
nature. 

This gap between what is known for the solutions of NSE (all of which are super-critical) and 
the sufficient conditions for regularity, is one of the manifestations of the important role scaling 
plays in the NSE well-posedness problem. 

There is a growing consensus that functional and harmonic analysis techniques alone would 
not be sufficient to break the regularity problem by obtaining a super-critical condition for well- 
posedness. Specifically, there are examples on NSE-like systems that blow-up in finite time despite 
many functional properties characteristic to NSE ([11, 18,27, 30]). 

This suggests a necessity of developing new approaches to understand NSE and non-linear 
systems in general. In particular, our results suggest that the stochastic multiplicative cascade 
framework introduced by Le Jan and Sznitman ([21]) may provide new insights into the NSE 
regularity problem. Indeed, in this note we establish a new scaling-critical condition for uniqueness 
of solutions, as well as provide evidence of a connection between the issue of uniqueness and 
natural scaling of the NSE. 

A natural way to explore the role of scaling in the theory of NSE is to consider scaling-invariant 
or self-similar solutions, i.e. the solutions satisfying 


u r = u , Pr = Pi Vr > 0. (1.3) 

Leray [25] observed that if u,p is a self-similar solution to (1.1), then upon choosing r = 
r(t ) = 1/Vi for fixed t > 0, one has 


u(x,t) 



X 

vi 


,i) 




where 

— AU — iu — ^(X ■ V)U + (U ■ V)U = —VP, V-U = 0. (1.4) 

Leray himself had the idea to use this self-similarity (backwards in time, with r(t) = l/y/T — t ) 
to produce an example of blow-up in the NSE problem. This was eventually proved impossi¬ 
ble due to the work of Tsai ([32]), as well as Necas, Ruzicka, and Sverak ([28]) (see also [24]) 
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who established that backward-in-time the only self-similar solution is 0. Study of forward in 
time self-similar solutions, particularly of equation (1.4), revealed several important existence and 
uniqueness as well as regularity results ([9, 14, 15,26]). In particular, Meyer ([26]) provided a 
framework of constructing solutions that are unique in a weak /."’-space starting from “small” ini¬ 
tial data. We note that the self-similar solutions must invariably posses singularity at the origin, 
as they are homogeneous functions of degree -1. The weak-L 3 space is a natural functional space 
for such functions. Later Grujic ([15]) showed that the solutions built by Meyer are in fact smooth 
(outside the origin). More recently, Jia and Sverak ([16,17]) proved existence of smooth solu¬ 
tions for (1.4), without a smallness assumption, pointing to a potential for lack of uniqueness of 
self-similar solutions for ‘large’ initial data and showing a pathway of how such solutions might 
be used to produce blow-up in Navier-Stokes equations. Cannone and Karch ([10]) also argued 
for the connections between the theory of self-similar solutions with large initial data and possible 
emergence of singularities in the NSE. 

The fact that self-similar solutions could be used to prove/disprove well-posedness for general 
NSE is another manifestation of the particular importance of scaling symmetries in the Navier- 
Stokes equations. 

1.2 The question of symmetry breaking - the description of the main results. 

In this paper we seek to provide an approach to both self-similar, as well as general NSE problems 
that could shed light into the specific issue related to this natural scaling, to be referred to as 
symmetry breaking, namely: 

Is the uniqueness of solutions to NSE tied to the uniqueness of self-similar solutions? 

If the solutions to (1.4) are unique, yet the solutions to (1.1) are not, then we have a manifes¬ 
tation of symmetry breaking in NSE, signaling that a possible lack of well-posedness could be the 
result of a mechanism that magnifies/creates deviations from natural scaling present in the initial 
data. On the other hand, if the uniqueness for (1.4) is closely tied to the uniqueness in (1.1), then 
the well-posedness problem is essentially connected to the natural symmetries of Navier-Stokes 
equations. Thus, the notions of scaling invariance and self-similarity, considered from the perspec¬ 
tive of symmetry breaking, provide the central focus of the present paper. 

We consider this issue in the framework of Le Jan-Sznitman stochastic multiplicative cascades, 
developed in [21], combined with the idea of majorizing kernels, introduced in [3], to investigate 
existence and uniqueness of a mild solutions to the NSE (see (2.2) below). In this framework, a 
multiplicative cascade process is associated to the mild formulation of NSE in Fourier space, and 
a solution is recovered form the initial data via an expected value of a certain recursive product 
along a generated tree. The space of initial data allowed is in part governed by the choice of the 
majorizing kernel (see Section 2 for details). 

In order to guarantee finiteness of the tree, a thinning procedure is usually employed. The 
thinning, which involves a chance of artificially terminating a branch, is guaranteed to generate 
a finite cascade, producing a unique mild solution to NSE, but at an expense of shrinking the 
smallness condition on the initial data. 

In contrast to the classical Le Jan-Sznitman approach, we will not employ a thinning procedure 
to terminate the cascade (see Section 2). Elimination of thinning is a step towards accommodating 
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wider families of initial data by relaxing, and eventually removing the aforementioned smallness 
condition. However, in the absence of thinning, one has to deal with a possibility of the formation 
of infinite cascade trees in finite time - the phenomenon called explosion , and our main object of 
study. 

In particular, we will not be concerned with the issue of existence of mild solutions built with 
such procedure (the solution is guaranteed to exist as long as the cascade is non-exploding and the 
associated expected values are finite - see Section 2). Also, we will not study regularity properties 
of such solutions (a difficult question, especially for more general spaces of initial data). Instead, 
our goal is to show that the explosion phenomenon in such cascades can be used as a surrogate for 
uniqueness for the solutions of the NSE in a certain functional class, allowing us to classify the 
associated uniqueness problems by the corresponding explosion time random variables. 

Specifically, we use this approach to study two families of NSE initial data: one governed by 
the Bessel majorizing kernel (2.6) (in which case we are able prove the non-explosion), as well as 
the dilogarithmic kernel (2.5) (which allows for a much wider space of initial data, but with a more 
nuanced explosion problem). 

We also adapt this approach to the study of scaling-invariant solutions (see Section 3). This 
forces a very different choice of the scaling parameter r (see (3.1)) than the one in (1.3), partly 
because the problem is posed in the Fourier setting. Nevertheless we show that the self-similar mild 
formulation we use - (3.2), and the Leray equation (1.4) are in fact equivalent (Proposition 3.4). 
Moreover, although the resulting cascade is quite different from the general NSE case described 
in Section 2, the dilogarithmic density appears naturally in the context of scaling-invariance (see 
(3.2)). 

The explosion problems themselves are defined in terms of the explosion time random variables 
in both non-symmetric and scaling-invariant cases (see Definitions 2.1 and 3.1) - critical (with 
respect to the scaling) quantities. Using the Le Jan-Sznitman martingale argument, we show that 
in the case of general NSE, the non-explosion of the associated multiplicative cascade provides a 
scaling-critical sufficient condition for uniqueness - see Proposition 2.1 and Remark 2.2. 

A natural question is to compare the case of dilogarithmic majorizing kernel in general, non- 
symmetric setting to the scaling-invariant case. While we were unable to fully resolve the associ¬ 
ated explosion problems, the main conclusion of this analysis - see Theorem 3.1 - is that at the 
level of cascades, 

The explosion problem is the same in both self-similar and general case (in dilogarithmic settings). 

This result provides evidence for a lack of symmetry breaking in the Navier-Stokes problem. 

The rest of the paper is organized as follows. 

In Section 2 we will use the Le Jan-Sznitman cascade without thinning, together with the 
idea of majorizing kernels, to formulate an explosion problem (Definition 2.1) closely connected 
to the issue of uniqueness of solutions to the mild NSE formulation (2.2) - see Proposition 2.1. 
Two particular kernels are considered (see (2.6) and (2.5) ). In the case of Bessel kernel, hb, we 
can prove the non-explosion (see Theorem 5.1 in the appendix). In the case of less restrictive 
dilogarithmic kernel, ha, we prove that the explosion is related to the uniqueness property of a 
certain solution to a non-linear PDE (see Proposition 2.2). 

In Section 3, an analogous procedure will be employed to arrive to an explosion problem in the 
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self-similar case (Definition 3.1). In particular, we relate the solutions obtained with this method 
to the solutions of the above-mentioned Leray equation (Proposition 3.4), and prove (Theorem 
3.2) that the explosion itself is a zero-one event (i.e. it is essentially deterministic). We also show 
in Theorem 3.1 that this explosion event has the same distribution as the explosion in the diloga- 
rithmic case described in Section 2 - the evidence towards similarity between the two uniqueness 
problems. 

Section 4 is devoted to comparison of the general and self-similar cases from the point of view 
of the symmetry breaking question, and discusses some open problems. 

Finally, Section 5 is the Appendix containing the proofs of the technical results related to the 
Bessel and dilogarithmic random walks, which appear in the functional settings adopted in Section 
2. It is worth noting here that dilogarithmic random walk (which arises naturally in the context of 
this paper) appears to be a new multiplicative stochastic process that may be of broader interest, 
e.g, see [19,22,23] for other occurrences of the dilogarithmic distribution in physics. 


2 Navier-Stokes Cascades & An Explosion Problem 

Next, we will describe the mathematical framework we use to analyze the existence and uniqueness 
problem for the Navier-Stokes equations more precisely, by specifying the meaning of “solution”. 
Due to lack of existence results for smooth (classical) solutions, the notion of solution in the weak 
sense is frequently used, where derivatives are in the distributional sense, as this allows one to 
search among functions u that are locally square-integrable in space. Namely, a divergence-free 
vector field u(t) is called a weak solution if for all v : M 3 —> M 3 - smooth, divergence-free 
functions with compact support, 


t 

(u(t), v) - (u(0), v) = J (z/(u(s), Av) + (u(s), (u(s) ■ V) v)) ds , 

0 

under the implicit assumptions on u that make the integrals above valid. Here, (u, v) = f R3 rr ■ v c/x 
is the (complex) L 2 inner product. This definition was introduced by Leray to provide a mathemat¬ 
ical framework that would accommodate the possibility that velocities may not be smooth at some 
“small set’’of points where “turbulence” is present. Indeed, Leray’s approach is proven to produce 
solutions that are are smooth except possibly a singular set of one-dimensional Hausdorff measure 
zero (see [8]). 

Taking Fourier transform in x in the equation above, we notice that 

t 

(u(t),v) - (u(0), v) = j (-i/(|£| 2 u(s),v) + (u(s),.F[(u(s) • V)v])) ds, 

0 

where w(£) = J r [te](^) = (27 t)~ 3 / 2 J r3 w(x) e~ tx '^d^ is the the Fourier transform of w. 

Using the reality condition u(—£) = u(£), the last term can be written as: 

1 —X 

(u(s),.F[(u(s) • V)v]) = ) 3/2 (u{s),u k {s) * (i&v)> = /2 (£ k u k (s) * u(s),v), 
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where v * w(£) = f R3 v(£ — rj)w{rj) drj is the convolution of functions v and w. Moreover, to 
incorporate divergence-free property in the first term of the inner product above, we can write 
-i i£i r 

( ( 27 r ) 3 / 2 ^ fc ^) * Q ( g )’ V ) = ( ( 2 7 t ) 3 / 2 J 

R 3 

with 

v(t?i) w(r7a) = -i(e$ • w(r/ 2 ))7r^v(r/i), (2.1) 

where = £/|£| and vr^±v = v — (e^ ■ v)e^ is the projection of v onto the plane orthogonal to £. 

Since |£| f R3 u(r), s ) ©^ u(£ — 77, s) dj) is divergence-free, we conclude that weak solutions 
satisfy 


u(t) - u(0) = J | -z/|^|-u(s) + J^3 /2 J U (77, s) © € u(£ - 77, s) dr) j ds , 

0 \ R 3 / 

which leads to the following mild formulation of the Navier-Stokes equations: 

t 

u(£, t) = u(£, 0 )e~ v ^‘ t + J e~ v ^ a J 6(77, t - s) © € u(£ - r),t - s) dr) ds . (2.2) 

0 R 3 

We note that weak solutions automatically satisfy the above mild formulation, and the solutions 
to ( 2 . 2 ) are weak solutions provided they are (uniformly locally) square integrable in both space 
and time variables ([24]). 

The stochastic cascade framework in Fourier space was introduced in [21] for the analysis of 
(2.2). The basic ingredients of the recursively defined stochastic object (cascade) associated with 
the problem (2.2) consists of (i) a continuous time binary branching Markov process in three- 
dimensional Fourier wavenumber space, and (ii) an algebraic operation ©^ defined in (2.1). The 
stochastic process is initiated with a Fourier mode 0 7 ^ £ g M 3 , where it holds for an exponentially 
distributed length of time Tg with intensity v |£ | 2 • Upon expiration of time T$, the particle either 
dies or splits into a pair of frequencies (modes) (Wi,W 2 ) G K 3 xl 3 . The random events of either 
dying or splitting occur with equal probabilities and independently of T$. In the case of a split the 
new frequencies are subject to the local conservation of frequencies condition 

W 1 +W 2 = ^, (2.3) 

and distributed according to 

E/(W t , W 2 ) = f f(v, t-v) '^r%Ts ~ d V. (2-4) 

J R 3 d * a(4) 

where h : M 3 \{ 0 } —* (0, 00 ) is a measurable function with full support for which h * //(£) < 00 
for each £ 7 ^ 0 , introduced in [3] as a majorizing kernel for (1.1). The two special choices given 
in [ 21 ] involved 

h d (£) = ^ 0, (2.5) 

h b ($) = e-M/\Z\, ^0, (2.6) 
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for which one has the identities 

h*h(Z)=c\Z\h(Z), (2-7) 

with c = Ci = 7r 3 if h = h d , and c = c 2 = 2ir if h = h b . We will refer to h b as the Bessel 
majorizing kernel, and to h d as dilogarithmic majoring kernel. The connection with Bessel and 
dilogarithmic distributions will be given in Propositions 2.4 and 2.3 (see also Remark 2.6). 

We define the conditional densities: 

hdjyphdte - rj) = |g| 

hd*hd{£) tt 3 |C - ^| 2 |^| 2 ’ 
h b (rj)hb($, — rj) 

h b * h b (£) 27r \rj\\$, - rj\ 

If the particle dies, the (Fourier transformed) forcing term, evaluated at its parent mode and 
appropriately scaled, is attached to the terminal node. Otherwise, if branching occurs this rule is 
repeated from each of the nodes at respective frequencies W) and W 2 . The focus of the present 
article is the unforced (g = 0) equation (1.1), in which case such a convention may be viewed as 
a “thinning”operation, that may not be necessary so long as there are only finitely many branches 
by any finite time t. If thinning is applied, however, then the associated genealogical tree is that of 
a critical binary Galton-Watson process and therefore in fact almost surely finite, e.g., see [1,2]. 

To be clear, in the case of no-forcing (g = 0) the Le Jan-Sznitman algorithm results in a thin¬ 
ning of the full binary tree that may be ignored so long as the branching process is non-explosive. 
This observation will be elaborated upon as a point of focus in the present paper. 

The algebraic operation 0^ is applied to a vector-valued function of the offspring (Wi,W 2 ) G 
M 3 x M 3 , provided by the initial data u 0 : M 3 —* M 3 , at each node of the genealogical tree having 
parental wavenumber £ as defined by 

u 0 (»7i) ©$ u. 0 ( 772 ) = -*(e$ • uo(r72))7r 4 ±u 0 (r7i), (2.8) 

where = £/|£| and 7r^±v = v — (e$ ■ v)e^ is the projection of v onto the plane orthogonal to 
Figure 1 shows a geometric interpretation of (2.8). 

This stochastic cascade provides a weak solution to (1.1) for initial data u 0 as an expected value 
of a cascade product under the algebraic operation ©$, 

u(tt) = |C|- 2 E € X( n ), (2.9) 

where X(r*) refers to a ©^-product of initial data and forcing at wave numbers determined by 
the branching Markov chain over nodes of the genealogical tree r t at time t, provided that the 
indicated expectations exist. The latter existence of expected values is an essential proviso whether 
the cascade is thinned or not. 

Of course the indicated expected values are to be interpreted component-wise when applied to 
vector quantities. In addition to the obvious decay required on the magnitude of the algebraic mul¬ 
tiplications for existence of expectation integrals, the rate of growth of the tree is also a significant 
issue. Specifically, we will be interested in the possibility of an explosion event in which infinitely 
many branchings occur within finite time. 

The phenomena of “explosion”of Markov processes and its relationship to uniqueness/non¬ 
uniqueness of solutions to the corresponding Kolmogorov equations is well-known in the theory 


Hi(<l 10 = 

Hdn I f) = 
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U 0 (*72) 


(e* ■ Uo(^2))7T 4 xu 0 (r/i) 


Figure 1: Geometric interpretation of uo( r 7i) 0$ no(r/ 2 ). 


of stochastic processes; e.g., see [2, 13,29]. The essence of explosion is that the stochastic process 
may leave the state space in finite but random time (, and then be instantaneously returned to 
the state space according to some arbitrary distribution. Moreover this regenerative extension can 
then be repeated to obtain a Markov process whose transition probabilities also satisfy the same 
Kolmogorov (backward) equations. However, this is a linear Markov process theory that does 
not directly apply to (1.1). Nonetheless, the associated branching process is a Markov process 
for which, in the absence of thinning, explosion cannot a priori be ruled out. In this context we 
consider the following. 

Definition 2.1. The explosion time of the Fourier mode cascade genealogy originating at £ 0 is the 
(possibly infinite) random variable given by 


n 



where for each n > 1, |s| = n, denotes a genealogical sequence s = (si, ...,s n ) G {1,2}", and 
for j < n, s\j = (s ],..., Sj) is the restriction of s to the first j generations. The random variables 
T s , s G U“ =1 {1, 2}", are i.i.d. mean one exponentially distributed random variables independent 
of the Fourier modes W s , s G U^ =1 {1, 2 } ”. The event [£ < oo] is referred to as an explosion event. 


Define also 


n 



and note that by monotonicity one has ( = („. 

Various general conditions for explosion/non-explosion will be given in the next section. A 
consequences of non-explosion is as follows. 

While the full connection is incomplete, the following provides some further evidence of a 
connection between explosion and the uniqueness problem for (1.1). 

Proposition 2.1. If there is no explosion then the stochastic cascade solution provides the unique 
mild solution to Navier-Stokes equations whenever the indicated expectations exist. 
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Proof. If explosion does not occur then the stochastic cascade is recursively well-defined and the 
same martingale (inductive) arguments of [ 21 ] may be applied. □ 

Remark 2.1. For the converse case that explosion can be shown to occur, one may construct 
the regenerative extensions as mentioned above. However it is not clear whether or not this is a 
pathway to non-uniqueness. 

Remark 2.2. Note that if we re-scale NSE according to seeding (1.2), then the explosion time 
random variable for the re-scaled cascade, ( r , and (from Definition 2.1 have the same distribution. 
Thus the non-explosion provides a critical (or seeding invariant) condition for uniqueness of the 
Navier-Stokes equations. 

From the point of view of uniqueness of solutions to associated PDEs, it is interesting to note 
that the absence of explosion does correspond to the uniqueness of solutions for an evolution 
equation. In order to state this equation we first define an operator A by 

A(/)(x)=^- 1 (|^|/(0)(x), (2.10) 

where J~ 1 denotes inverse Fourier transform. 

This operator acts to increase the higher frequency oscillations of /(x) by the same magnitude 
as differentiation. As this operator is closely related to differentiation, it is known as a pseudo¬ 
differential operator. The evolution equation associated to the branching process as follows, which 
contains this pseudo-differential operator, is as follows. 

Proposition 2.2. Let h = h^or h = hi, and assume that the pseudo-differential equation 

f ^ Ar • cN(v 2 ) 

\ v(x, 0 ) = .F _ 1 (/i)(x) 

with c — n 3 or c = 2" respectively, has a unique mild solution satisfying i;(£. t)\ f //(£) for edl 
t ^ 0. Then explosion does not occur, i.e. P([£ < oo]) = 0. 

Proof. Fix h = h,i- or h = hb and c = tt 3 or c = 2 tt respectively. Let Z((, t ) denote the number of 
offspring by time t > 0 for the Fourier mode cascade genealogy starting at time t — 0 at frequency 
0^6 M fc . In particular, after an exponentially distributed time with parameter \(\ 2 . the parent 
particle is replaced by two particles of frequencies 77 , ^ — 77 where 77 has probability density 

H(ri\£) = h(t-Ti)/(ch*h(t)). ( 2 . 11 ) 

For k ^ 1, let 

m((,t,k) = P^Z(tt)=k). (2.12) 

In particular m(£, t, 1) = exp(—|^| 2 f). Explosion does not occur if and only if P^(Z(^, t) = 00 ) = 
0 , or equivalently YlkLo rn ((i T k) = 1 . 

Let v(£, t , k ) = /i(^)m(^, t , k ). Then trivially v(fi, t , 1) = h (£) exp(— \$,\ 2 t). Moreover, condi¬ 
tioning on the time of the first branching, it follows using (2.7) that for k ^ 2, v((, t, k ) satisfies 
the integral equation 

k—1 n. 

v(€,t, k) — c|£| y2 / / v(r],t-s,j)v(^-rj,t-s,k-j)drje^ 2s ds, v(£,0,k) — 0. (2.13) 
j=1 Jo Jr 3 
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Let v(£,t) = Y^k=i v (£-f k). This series converges since all terms are non-negative and the 
partial sums are clearly bounded above by h(£). Note also that 


H£,t) = h{g)^2m{£,t,k). (2.14) 

k =0 


Summing (2.13) for k ^ 2, and adding the missing term v(£,t, 1) one finds that v satisfies 

v(£,t) — h(£) ex;p(—|£| 2 f) + c|£| I f v(rj,t — s)v($, — rj,t — s)drje^ 2s ds (2.15) 

Jo Jr 3 

It follows that d(x, t), the inverse Fourier transform of v(£,t), satisfies the following reaction- 
diffusion equation of [27], 

dv 

— (x,f) = At>(x,t) + cA(u 2 )(x,t) (2.16) 

with initial data u(x, 0) = J r_1 (/r)(x). 

One may easily check that with h = 7 r ~ 1 (/r), 

Ah + cA(/i) 2 (x) = 0. (2.17) 

By hypothesis, u(x, t) = h(x) and thus, using (2.14), 


Pd C >t) = PdZ(S,t) < oo) = Y,™(S,t,k) = 1, V£,t. (2.18) 

k =1 

□ 


Remark 2.3. Note that since h,i and h b are radially symmetric, the cascade, hence m(£,t,k) 
and rn(f, t ) = fffk=o rn (£i T k) are also radially symmetric. In particular, it follows that with H 
defined by (2.11), 

m(|£|,t) = exp(— \$,\ 2 t) + |£| 2 f [ m(\ri\,t — s)m(\$, — r]\,t — s)H(rj\£) drje~^ 2s ds (2.19) 

Jo Jr 3 

Remark 2.4. In 2007 Chris Orum has announced the equation of Proposition 2.2, and its role 
in the explosion problem in a session of the 32nd Conference on Stochastic Processes and their 
Applications, Champaign-Urbana. However the uniqueness/explosion problems remain unsolved 
for general majorizing kernels h as initial data. 


We conclude this section with a small further elaboration on the probabilistic significance of 
the two kernels hd, hb- Additional features are discussed in the appendix. 

Let 

nm? „, m 

a(?) = hzwr ( ’ 


Then, for either kernel h = ha or h = hb, one has that af$,) = c _1 |£|, i — 1,2 (2.20) with 
c = 7T 3 or c = 27 t respectively, defines the same pseudo-differential operator given by a positive 
multiple of sj— A. However the following two propositions dramatically distinguish the associated 
branching Markov chains. 
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Proposition 2.3. Assume that h(£) = h d (£) = |£| 2 ,0/^e M 3 . Then for each s E {1, 2}°°, the 
sequence {|W s |j + i|/|W s |j| : j — 0,1,... }, W s | 0 = £, is an i.i.d. sequence under P,\, such that, 

IW„ 


Pa 


s|j+i| 


|W 


sb'l 


E dr) = 2n 2 In 


1 + r 


1 — r 


dr 

—, r > 0. 
r 


Pd ^ 


|W 


«b+ i i 


iw 


s\j\ 


E dt) = 2n 2 In | coth(t/2)| dt , t E 


Proof. Part (ii) is an immediate consequence of (i) by a change of variables formula. As noted 
above, the distribution of |W| was computed as (1.22) of [21]. Essentially the same calculations 
apply to the ratios of magnitudes as follows: For non-negative and integrable g on (0, oo), one has, 
since h * h(£) = vr 3 /|^|, and using (2.4) that 



71 3 |£l / 9 ( 777 


\m 


dr] 


7r 


-3 


£\J |r7| 2 |^ - 77| 2 

s(M)ot dV 


viu — V 


7r 


-3 


iKM)' 


2 ’ 

dv 


u = 


7r 


-3 


'0 J |w|=l 

/»00 /»7T 


2t r 


-2 


>0 JO 

poo 


2tt 


-2 


g(r) 
g(r) In 


v | 2 (|v| 2 — 2 u • v + 1) 
dw dr 

q(r)— - 

7 r 2 2r u • w + 1 

sin f d(f> dr 


r 2 — 2r cos 0 + 1 
1 + r| dr 


1 — r\ r 


□ 


Definition 2.2. The multiplicative random walk (|W S | J+1 | : j = 0,1,...} on (0, oo) will be 
referred to as the dilogarithmic random walk, or dilog random walk for short. 


Proposition 2.4. For arbitrary £ E M 3 , let W denote the random vector in R 3 with density 

e l£l e -H e ~\£,-n\ 


n b (r, 10 = 


2tt \rj\\£ — rj\ 


Then 


P(|W| e dr) = 


l^j-e 2r (e 2 ^ — 1) for r ^ |£| 


lid 


,-2 r\ 


for 0 ^ r ^ |£| 


In particular, along any path s, the sequence |W@| = |£|, |"W S |! |, |W S | 2 |, ..., is a Markov chain 
with stationary transition probability density 


P(u,v) = 


I e 2v^ £ 2u _ -Q f or v ^ u , U > 0 
-(1 — e~ 2v ) for 0 ^ v ^ u, u > 0. 
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Moreover, 


E*|Wi| = 


ICI + 1 


( 2 . 21 ) 


Proof. For arbitrary £ G M 3 \ {0}, let W denote the random vector in M 3 with density 

e l£l e -\v\ e -\€-v\ 


Ufa 10 = 


2 ^t |ry||^ — 17 1 


Let u = |£|, and use spherical coordinates with p = 1 77 1. Direct calculations exploiting the cylin¬ 
drical symmetry of the distribution give 


P(|W| > r) = —(27 t) / e~ p p 


r exp(-(p 2 - 2pucos(j) + u 2 ) 1/2 ) . 

f -/ 2 — o - , , 9 \ 1 /9 -sm fidfidp, 

0 (p z — 2pu cos 0 + u-) 1/z 


= -(-!) ^(e 


- p( P -(p+ u ) _ P -|p-“h 


u 

e u 

u 


dp, 


-pp-\p-u\ 


e 


dp -e 


-2 r 


For r ^ u, the integral in the last expression is (e“/2)e 2r so that in this case, 




P(|W|>r) = — (e 2 “-l). 


Similarly, for r ^ u one obtains that 


f*U /»oo' 

+ / 


e -p e -\p-u\ _ p -u 


dp = e U (u-r +-), 


so that in this case, 
and the result follow by differentiation. 


I II — V 

P(|W| > r) = —(1 — e~ 2r ) +- 

2u u 


□ 


Remark 2.5. The calculation of the marginal distribution of \ W, can be found in ([21], Proposi¬ 
tion 2.1) for h = h f ]. A similar calculation was provided here for ease of reference. The Markov 
property follows by exploiting the construction together with the form of the transition probabilities 
as functions of the norms; see ([2], pp 502-503). 


Remark 2.6. Euler’s dilogarithmic function may be defined by 

or j 

Li 2 (r) = — ln(l — u) — , r < 1. (2.22) 

Jo u 

The dilogarithmic function is a special case of polylogcirthmic functions Li s (.x) whose domain of 
definition may be extended to include complex values of both x and s. An extensive literature is 
available for properties and relationships between polylogarithmic functions, with connections to 
Bose-Einstein and Eermi-Dirac statistics, e.g., see [19,22,23] 


The explosion problem is solved for the Bessel kernel in the appendix (Theorem 5.1), but it 
remains quite illusive for the dilogarthmic kernel. However, as will be seen in the next section, the 
dilogarithmic kernel is somewhat singled out by the self-similarity cascade. 
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3 Self-Similar (Navier-Stokes) Cascade & And Its Associated 
Explosion Problem 


In this section, we obtain a stochastic cascade associated to the Navier-Stokes equation when self 
similar solutions are considered. It should be remarked from the outset that the kernel Hd occurs 
naturally in this situation, as a direct consequence of the scaling properties of the solutions of the 
Navier-Stokes equations. We present first the mild formulation of the Fourier transform of the 
Navier-Stokes equations for self similar solutions. A probabilistic representation for the solution 
of the resulting equation is given in terms of what we call the self similar cascade. We show that 
important statistical properties of this self similar cascade and the Navier-Stokes cascade obtained 
using the dilogarthmic kernel Hd are identical. 

As noted in the introduction, the scaling invariance of the Navier-Stokes equations show that 
if u(x, t),p(x, t ) is a solution then for any r > 0, u r (x, t) = ru(rx, r 2 f), p r (x, t) = r 2 p(rx, r 2 t ) 
is also a solution of the Navier-Stokes equations. Assuming the initial data is also scale invariant, 
uniqueness would imply the self-similarity property u r = u. In the Fourier domain, this scale 
invariance corresponds to, with v = u r , 

v(&t) = W-yt), 

ry* — /vi 

so, with r — |£|, 

v(£,f) = |£| 2 £) ( 3 . 1 ) 

where = £/|£|. Thus, since v satisfies (2.2) it follows that a self similar solution of the Navier- 
Stokes equations satisfies 

u(e 6 |£| 2 t)) = e“ t|c|2 Uo(e € ) 

+ (2tt)“ 3/2 [ e _l4|2( *“ s) |^| 3 [ u(e r? , |? 7 | 2 s) u(e € _„, |£ - T 7 | 2 s) 1 - -d'qds. 

Jo J |»7| Is — ^1 

The change of variables r] = l^ry', s' = |£| 2 s, and with A = \$,\ 2 t, gives 


u(e s , A) = e u.jic.j 
+ (21)- 


CO 

1 

•< 

T 

03 

K __ 

f u(e,j, \rj\ 

Jo J 



\rj\ z \e£ — -q | 


drjds. 


Recall that Hd(rj |e^) = ( 7 r 3 |? 7 | 2 |e^ — 77 I 2 ) 1 so one has 


e (A s) / u(e^, |r 7 | 2 s) 0 ^u(e e? „ T? , \e^-r]\ 2 s)H d (r]\e i )dr]ds. 

(3.2) 

We refer to the parameter A > 0 as the similarity horizon. 

A probabilistic interpretation for (3.2) follows similar steps as those introduced before. Con¬ 
sider a binary tree rooted at 0 with vertices indexed by V = U n ^i{l, 2} 11 - see Figure 2 for an 
illustration. Denote by dV = {1,2} N . Elements in each of these sets are denoted by s and (s) 


u(e^,A) = e A u 0 (e € ) + (7r/2) 3/2 [ 

Jo 
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respectively. Let {T s , s G V} be a collection of i.i.d. random variables with an exponential dis¬ 
tribution with parameter 1. Given a direction e s , let W s i be a random variable with distribution 
Hd(rj\e s ), set W s2 = e s — W. s i and for j = 1, 2, define the directions e S j = W, S j/1 W. SJ |. Finally, 
given a horizon X s , define for j = 1,2 A S j = |W SJ | 2 (A S — T s ). On each (s) G dV , the branching 
process stops at level 

N {s) = inf{m ^ 0 : A (s | m) < T {s \m)}- (3.3) 

Completely analogous to (2.9), the solution of (3.2) is then given as an expected value of a 
recursive product involving the algebraic operation Q e$ provided this expectation is finite. Fur¬ 
thermore, the evaluation of this recursive product can be done if and only if along any path in the 
binary tree, the random variable N( s ) defined in (3.3) is finite. 

From the definition of the random variables one has 


A( s |n) - T {s | n ) = ((...(((A 0 - T 0 )|W( S | 1 )|- - r (s |i ) )|W (s | 2 }| 2 - T( s | 2 ))...)|W (s | n) | 2 - T {s | n) ) 


(niw*i 2 )L.-E^> 

k =0 V j =0 


n j 


k =0 


IW 


(s\k) | 


where we have used that |W( S | 0 )| 2 = 1. Thus, for given A 0 and ( s) G V, the event [N^ 
equals the event 


inf{m ^ 0 : 

3=0 11 

This motivates the following definition. 


l=o 


W 


^ A 0 } = n. 


(s\k )| 


n 


Definition 3.1. For a fixed unit vector e 0 , the similarity explosion horizon is the (possibly infinite) 
random variable 


CM 


lim inf > T/ ,i,-\ 

n —>oo |s|=n 


3 

k =0 


w 


(s\j) I 


The self similar explosion event is defined as A eo = U m j>i [C( c o) < rn ] so that PL4 C(I ) is the 
probability of self similar explosion. 


Note that with 


CM) = inf 
\s\=n 


n 


/ . T{ s \j): 

3=0 


m =0 iw 


<s|i> I 


one has, by monotone convergence, that C( e o) = lim n _>.oo Cn(eo)- 

While the self-similar cascade construction is quite distinct from that of the Navier-Stokes 
cascade, one may note that for fixed (s) G <9V, the random variables Rj = \ W < a |_,-)\,j ^ 1 are i.i.d. 
with the dilogarithmic distribution with density 


V[r) 


7T 2 r VI 1 - r l/ 


Indeed, since the distribution of W( s m depends only in the unit vector e 0 , the proof of Proposi¬ 
tion 2.3 shows that Ri has the dilogarithmic distribution. The claim for R 3 follows by induction. 

In order to relate the explosion problems for the self similar cascade and the Navier-Stokes 
cascade, we have the following result. 
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Proposition 3.1. For any n f 0 , the distribution of Q, is independent of the initial direction and 


Cn( e o) inf VI T( s \j) „ 

\ s \= n ~~Z J 1 


3 =0 


J 

k =0 


Rj I 


where Rq = 1, and {/? 7 , A a sequence ofi.i.d. random variables with density V{r). 

Proof Let Q be an orthogonal 3 by 3 matrix and e a unit vector in M 3 . Let r). rf be random vectors 
distributed according to Hd(rj\Qe) and //,/(r/|c) respectively. It follows easily that in distribution, 
f/ and rf are equal and thus independent of the particular initial direction e used in //. The proof 
is completed, since as noted above, V{r) is the density of If. □ 

As a consequence of Proposition 3.1, the distribution of the sequence \{ v \j),j ^ 1 is also 
independent of the initial direction e 0 . 





A 2 2 < T 2 2 


Figure 2: Self-similar cascade with explosion cartoon. 


Moreover, comparing the Q, above with ( n - defined in the context of Definition 2.1 for the 
kernel h ( i - we obtain our main result connecting the self-similar and dilogarithmic uniqueness 
problems. 
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Theorem 3.1. The events [Cn(|£|) > t\for a dilogarithmic density and [Cn( e ^) > t\€\ 2 ] have the 
same distribution independent on the choice of the initial wavenumber £ or initial direction e^, 
and hence the explosion time ( from Definition 2.1 for the dilogarithmic kernel, and the eventual 
explosion (from Definition 3.1 have the same distribution, independent of the choice of or e 0 . 


Proof. Recall that when the dilogarithmic kernel is used to determine the distribution of the 
branching frequencies, it follows that for any ( s) € dV 


Rk 


|W(.| fc> | 

W (s | fc —i) 


k^l 


is a sequence of iid random variables with density V{r). Now, with W s | 0 = one has 


C«(l*l) = 


v 



inf T, 


M= n 


s| 0 


En 

j =1 k =1 


Wjjt-iP 

|W s|l |2 * b 




1 


n 


3 

k =0 


R k I 2 


where R.q — 1. Thus the we obtain the equality, in distribution, of the events \Q n (|£ |) > I ]. and 

[Cn(e € ) > tm i 

In analogy with Proposition 2.10, one has the following 


Proposition 3.2. Let 

OO 

Z{ Ao) = l + ^^i[T 1 , <AJ. 

n=0 \v\=n 

Define m(X, k) = P(Z(A) = k).Letrh( A) = YlkLi di(X, k). Then, 

p\ P7T POO 

m(A) = e _A + / / / m(r 2 s)m((l — 2r cos 9 + r 2 )s)H(8, r) drdOds. (3.4) 

Jo Jo Jo 

Moreover, if r'n(X) = 1 is the unique non-negative solution then there is no similarity explosion. 

Proof. Note that Z(X 0 ) represents the number of branches of the self-similar branching process 
started with horizon A 0 . Recall that from the definitions of A iS , Z is independent of the initial 
direction. Each time the indicator does not vanish, a branching occurs increasing the number of 
branches by 1. The extra term is to count the initial branch. 

For k f 2, condition on the time of the first branching to get, 



m(\r)\ 2 s, j)m(\e 0 - r/\ 2 s, k — j)H d (r/\e 0 ) dr)ds, 


where e 0 is arbitrary. 


(3.5) 
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Let H(9, r) denote the average of H d with respect to the angle 0, 

H(0 r) = 2 Shld 

{ ’ f 7T 2 1 — 2r cos 9 + r 2 ' 

Then, the independence of the above equation with respect to the direction e 0 is better illustrated 
in the following; 




m(r 2 s,j)m(( 1 — 2 r cos 9 + r 2 )s , k — j)H(9, r) drdOds. 


Summing on k the previous equation and adding the term corresponding to k = 1, one has (3.4). 
It is clear that rh = 1 is a solution of this equation, so non explosion is equivalent to showing that 
this is the only non negative solution that is bounded by 1. □ 


While we can not prove that m = 1 is the only solution of (3.4), we note that the behavior at 
infinity can be used to determined if fh{ A) < 1 on a set of positive measure. In fact, if for some 
e > 0, ?n(A) ^ (1 — e) on a set E of positive measure, then fh is bounded by a decreasing function. 
Indeed, for any A > 0, 0 ^ rh( A) ^ 1, and from (3.4) one has 


m(A) ^ e A + / e 


3 -(A-s) 


m(r 2 s)V(r) dr 


< e A + / e s ds — e e s V{r) drds 


= 1 — e/i(f?)( 1 — e A ) 


where fi(E) = J E V(r) dr. 

We are now ready to establish one of the main results of the paper. Define the finite horizon 
probability of explosion in a similar way as that of self similar explosion. To be precise, let 


A(A) — n n>1 [Cn ^ A]. 

Then, rh( A) = 1 — P(A(A)), and thus 

limsupm(A) = a < 1 P(z4(A)) > 0. 

A—^oo 

We then have the following; 

Theorem 3.2. The self similar explosion event is a 0,1 event and independent of the initial direc¬ 
tion. 

Proof. Note that since P(A(A)) is an increasing function of A, m(A) is decreasing. Let 0 ^ o ^ 1 
be defined by liniA^oo fh{\) = a. Using dominated convergence, one can take limit as A —* oo in 
(3.4) to get a = a 2 , so a = 0 or 1. The proof is completed, since P(/l eo ) = 1 — a independent of 
e 0 . □ 


An important consequence of this result is the following corollary. 
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Corollary 3.1. For any £ f 0, the explosion event for the Navier-Stokes cascades defined in terms 
of the dilogarithmic kernel Fid is a 0, 1 event. 

Proof. The corollary follows from the equality, in distribution, of the events [Cn(|£|) > t] and 
[C(e^) > t \£\ 2 ] (See Theorem 3.1). □ 

Similarly, the integral equations (2.19) and (3.4) can be shown to be equivalent in the case the 
Navier-Stokes cascade is defined using the dilogarithmic distribution. 

Proposition 3.3. Let m(|£|, t) be the solution of the integral equation 

m(|£|, t) — e - ^ 2 * + |£| 2 [ f m(\r]\,s)m(\£ — rj\,s)H d (ri\£) drjds. (3.6) 

Jo Jr 3 

Then 

rh(X) — m(|£|, A/|£| 2 ) (3.7) 

is a solution of (3.4) Conversely, given a solution fh(X) of (3.4), equation (3.7) defines a solution 
of (3.6). 

Proof. Introduce new variables rj = = |^| 2 s, and recall that H d (rj\£) dr\ = H d (rj'\e^) dr/'. 

Then changing variables in (3.6), and with A = |^| 2 t, one has, 

A/|C| 2 ) = e" A + [ e“ (W) [ s'/\^\ 2 )m(\^\\e^-rj'\, s'/\^\ 2 )H d (rj'\e^) dri'ds'. 

Jo Jr 3 

With rh( A) as defined in (3.7), one has, dropping primes, 

m(A) = e _A + / / rh(s\rj\ J )m(s|e^ — r/\ 2 )H d (ri\e^) drjds. 

Jo Jr 3 

The proof is completed by noting that (3.4) is obtained from this equation by integrating the angular 
variables and, to obtain the converse, reversing the steps. □ 

3.1 Self Similar cascades and Leray equation 

In this subsection we show that the self similar stochastic cascade can be obtained directly from 
the Leray forward equations (1.4) 

Proposition 3.4. Let U(X) be a solution of the Leray equation (1.4), U denote its Fourier trans¬ 
form. Then, with a unit vector in M 3 and A > 0, 

u(e$,A) = AU(A/Ae^). 

satisfies (3.2). In particular u(e^, 0) = uo( e ^) 


19 


Proof. Recall that the forward Leray equations are obtained assuming a solution of the Navier- 
Stokes equations of the form 


u(x,t) = y^U(x/v / f), 


and are given by 


-AU-^U- i(X ■ V)U + (U • V)U = -VP, V • U = 0. 
Taking Fourier transform and projecting on divergence free vector fields, one gets 


(1 + |£| 2 )U+^- V)U + (2tt)- 3 / 2 | 


U(£ — ry) ©£ U(r 7 ) drj = 0 


(3.8) 


Let = £/|£| and define V(e^,r) = U(re^). Since 


«• vu = rd f, 

dr 


one has, with some abuse of notation, 


(1 + r 2 )V + -r^- + (27r) 3 ^ 2 r 
2 dr 


U (£ - r ]) U(77) dr ] = 0. 


Multiplying the equation by 2re r , one obtains 

-f-frW) = -(2 tt)- 3/2 2 rV 

dr 


U(^ - 77) ©£ 11(77) dr ]. 


Let V(e, r) = r 2 V(e, r). Then 

A( e Pv) = —( 2 tt )- 3/2 2 re r 
dr 


V (e re -r,, |re - 77 ]) 0 € V(e ?7 , |t 7 |) t - 

17 ™e — 77 1 ^ 1 77 1 


drj. (3.9) 


Note that one factor of r is used to get, up to a constant, H d (r]\re). 
One may easily check that 

lim V (e, r) = u 0 (e). 


r —>0 


Indeed, since u (t, £) = tXJ( \Ztf ), for $, = we have: 

Uo(et) = limu(Lee) = lim t U ( y/t , £) = limtVYe*, Vt) = lim V(e^, Vt). 
v s t->o w t->o t-s-o 5 t->0 * 

Integrating equation (3.9), and accounting for the constant to get H d , we obtain 


e r V(e,r) = u 0 (e) - (vr/2) 


3/2 


2 se s 


|se - 77I) ©^ V(e v , \rj\)H d (rj\se) dijds. 


With the change of variables 77 = sr/', and noting that e se _ v = e e - T y and that ^(77! se) dr] = 
H d (r]' |e) drj', we have, dropping primes 


e r 2 V(e,r) = u 0 (e) - (tt/ 2) 3/2 / 2se s2 / V(e e _^, s|e - 77I) ©* V^, 5177I ) iP rf (77 |e) dr] ds. 
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Let t = s 2 to get 

V(e,r) = e“ r 2 u 0 (e)-( 7 r/ 2) 3/2 [ e~ {r2 ~ t) I V(e e _ v , v^e-r/DQ^V^, Vt\rj\)H d (r]\e) dr] dt. 

Jo Jr 3 

The proof is completed setting A = r 2 and defining u(e, r 2 ) = V(e, r). 

□ 

Remark 3.1. As an aside, one may note that the choice of the scaling parameter r is completely 
arbitrary. Corresponding to the choices r = Ij\ft made by Leray, and say, r = l/|x|, respectively, 
let U!(x, t) = (l/v / f)U(x/v / t), and u 2 (a:, t) = (l/|x|)V(x/|x|, t/|x| 2 ) Let’s note that U and V 
can be related by an application of the Kelvin transform T\ with respect to the unit sphere in M 3 . 
To see this, recall that 7^[u(y)] = (a/|y|)u((a 2 /|y| J )y), defines the Kelvin transform of u with 
respects to the sphere of radius a. Now, letting X = x/s/i one has 

Ti[(iM)u(x)] = -A- u(x/|x| ! ) = 

-!-U(x\/i/|x| 2 ) = Tu(x/|x|, ^t/|x|) = Tv(x/|x|.i/|x| 2 ). 

|x| |x| |x| 

4 Conclusions and Further Directions 

The primary goal of this article was to precisely formulate a notion of symmetry breaking for 
the three-dimensional incompressible Navier-Stokes equations, and to provide an approach to the 
resulting symmetry breaking vs or not dichotomy. The notion that is introduced builds on a variant 
of classic scaling and self-similarity ideas of Leray [25]. Namely, symmetry breaking is defined 
as a phenomena in which one has uniqueness of self-similar solutions, but non-uniqueness of 
general solutions. The approach is derived from a stochastic cascade representation (NSC) of the 
Navier-Stokes equations introduced by Le Jan an Sznitman [21], together with a corresponding 
development of a cascade representation (SSC) of (mild) self-similar solutions. The essence of the 
approach is to exploit a notion of branching process explosion as a surrogate to non-uniqueness. 
A branching random walk cascade, namely the binary branching dilogarithmic random walk on 
(0, oo) viewed as a multiplicative group, is obtained as a common element of both representations 
for comparison. A principle result was the equivalence of the explosion phenomena for (NSC) and 
(SSC). In addition it is shown that the explosion criteria is critical in the sense of scaling, and a 
zero-one law is established for the explosion event. 

It remains to firm up the precise connection between explosion and non-uniqueness. A related 
semilinear pseudo-differential equation of Proposition 2.2 and an integral equation of Proposition 
3.2 can be associated with the branching numbers in such a way that uniqueness of solutions to 
either in an appropriate space is shown to be equivalent to non-explosion. In fact, although not 
obvious, as shown by Proposition 3.3, the two equations are equivalent. However the yet unproven 
connection between explosion criteria and uniqueness is expected to be that non-explosion corre¬ 
sponds to the uniqueness of mild solutions represented by (NSC) and (SSC), respectively. Proving 
this in appropriate function spaces is a substantial challenge to the overall approach. Assuming that 
this will be achievable, the surrogate results will prove that the equations are in fact not symmetry 
breaking. 
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5 Appendix: Bessel and Dilogarithmic Markov Chains & Ex¬ 
plosion 


This appendix records some general approaches to the explosion problem that may eventually 
prove useful as we learn more about the dilogarithmic branching random walk. In fact, we are able 
to demonstrate their effectiveness when applied to the simpler case of the Bessel kernel, which we 
show to be non-explosive. At a heuristic level, it is the mean reverting property (2.21) that makes 
the Bessel kernel simpler to analyze. 

The first approach to explosion exploits the monotonicity in the sequence {Cn}- 
Proposition 5.1. Let ( be as in Definition 2.1, and assume that for some A > 0, 

n A 

2n]E i€i n A+ i W .i2 °> as n °°- 

j =i 3 

Then P([C = oo]) = 0. 

Proof. To prove non-explosion it suffices to show that for any B > 0, 


or equivalently, that 


P\£\ (Cn > B eventually) = 1, 


flsidC < B]) = P| 5 |(n“ = 1 [C„ < B]) = lim /^(C„ < B) = 0 


where we have used the monotonicity of the sequence {Cn} - For the latter observe that, 


^41 (Cn < B) 


n 

.7=1 


< ^p m (Y.m b \- 2 T lu < b) 

3 = 1 

= 2 n P\e\(e- X Z?=i\W 1 \ i $*T llj > e -AB) 


for any A > 0, where 1 \j = (1,1,..., 1) is on the fixed, but otherwise arbitrary, tree path (1,1,...). 
By the Markov inequality, one has 


P m (C n <B) < 2 n %e -AE ?=i |wri w | 2 B\j e XB 


= 2v« E|£ ,n 


A 


l\\+\ Wlb r 


which converges to 0 as n —> oo. 


□ 


As an illustration of this methodology we provide a proof of Theorem 5.1 below, establishing 
that the branching Markov Chain defined using the Bessel kernel hb(£) to determine the distribu¬ 
tion of the branching Fourier frequencies does not explode. The mean reversion property (2.21) 
provides some indication as to why one may expect the corresponding branching Markov chain to 
be non-explosive, as will be shown is indeed the case. Namely, 
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Theorem 5.1. The explosion horizon is almost surely infinite for the Bessel Markov chain. 

For the proof we first note the following more refined property of the Bessel Markov chain. 
Lemma 5.1. Assume that W is a non negative random variable with probability density 

, , f -e~ 2w (e 2u — l) for w^u,u> 0 

Pu{w) - | l( 1 _ e -2«,) for 0^W^W,W>0. 


where u is an arbitrary positive constant. Then 

A 


E„ 


-A + 1L 2 

Proof. Use integration by parts to note that 
f A 


< 7ta/A, Vw, A > 0. 


A + w 2 


e 2w dw = y/~\ arcta.n(w/\/A)e 2w + 2-\/A / arctan(w/-\/A)e 2w dw. 


One has 

E ''' 


A 


' A + w 2 J o A + w 2 
1 f u A 


U In A + W‘ 


p u {w) dw 


dw- 

u 


- --e 2w dw + - 

A + w z u 


A 


A + w 2 


g— 2(id—m) dw 


1 1 

—— \/Aarctan('u/\/A)- 2y/\ / arctan(w/-s/A)e~ 2 "' dw 

U U In 


1 1 

— V^\arctan(w/\/A) H— 2s/~\ / arctan(w/\/A)e~ 2< ' 1 '' _n ) dw 


u 
2a/A 


u 


< 


u 

7 Ta/A 

u 


arctan(w/\/A)(e 2 ^ w — e 2w ) dw 


2s/\ 


u 


arctan(w/\/A)e 2w dw 


( e - 2 (—) _ e ~ 2w ) dw = ny/X — (1 - e- 2 “). 

2u 


The result follows by noting that (1 — e x ) / x f 1 for any x. 


□ 


Proof of Theorem 5.1 

Note that successive use of conditional expectations on the sigma field generated by the 
branching process up to the j th branching event, and Lemma 5.1 one has 


%n 


A 


j= 1 


A + |W;| 2 


= E 


1*1 


n— 1 


n A+ iwd 2 E ' w ^A + iwj 2 < 


u=i 


^ (7r\/A) E| C | 

^ {jty[X) n 


n —1 

n 

.3 = 1 


A 


A + IW, 


The theorem follows applying Proposition 5.1 with A < l/(27r) 


□ 
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Use of the monotonicity approach is less transparent for analysis of the dilogarithmic explosion 
problem. Another approach is generally possible that builds on a variant of the Biggins-Kingman- 
Hammersley (BKH), e.g., see [4-7], computation of the speed of the leftmost particle for additive 
branching random walks in terms of multiplicative branching random walk. It is potentially appli¬ 
cable to the dilogarthmic kernel precisely because for any path s G {1, 2}°°, 

, , , A |W sh -| 

W *|n = ffl lW l ’ W = 1 ’ 2 '-’ (5-1) 

f = \ l W S |j-ll 


and the ratios are i.i.d. That is, for any path s G {1, 2}°°, the sequence {|W s |j| : j — 0,1,... } is 
a random walk on the multiplicative group (0, oo) starting at |W S | 0 | = |£|. That is, for the diloga¬ 
rithmic kernel the branching Markov chain is in fact a branching random walk on the multiplicative 
group (0, oo). 

First let us recall the general heuristic underlying (BKH) speed calculations on the additive 
group of real numbers: Suppose that {S'™ : n = 0,1, 2,... } is an additive random walk on K with 
mean zero and starting at zero. Then by the weak law of large numbers P(S n < nc ) —>■ 0 as 
n —> oo for any c < 0. Let rri(0) = Ec ns ' and consider the following large deviation inequality 

m(d) n = Ee eSn 

> e nGc P(S n > nc). (5.2) 

Thus, 

P(S n > nc) < exp {—n(9c — \nm(9))}, (5.3) 

and in particular, 

P(S n > nc) < exp{— nsup(9c — ln?n( 6 1 ))} = e~ nI ^ c \ (5.4) 

c<0 


where 1(c) = sup c<0 (6 ) c — In rn(0)) is the Legendre transform of In m(6) at c. The Cramer- 
Chemoff theorem provides general conditions for which 


\\\P(S n >nc) 
hm - 

n—kx) n 


-1(c). 


To apply this to the computation of the speed of the left-most particle of a branching random 
walk one reasons as follows: At the n-th generation the expected number of particles located to 
the left of c < 0 is 2 n e ~ nI ^ for large n. Thus the extremal speed is given by 7 = c such that 
2 n e -ni(c) ^ p The (BKH) theorem confirms this. For the calculations involved here it is actually 
enough to calculate a lower bound on the speed. 

This principle translates to the multiplicative group as follows: 


Proposition 5.2. Consider a binary branching random walk on the multiplicative group (0, 00 ). 
That is, the n-th generation particle position for the genealogy s = (. 7 ...., s n ) G {1, 2} n is given 
by the product n”=i Y s \j, where (Ly*i, U* 2 ) ’s are i.i.d random vectors with positive components. 
Then 

n 

lim min(TT T' s i 7 ) 1//ri = e 7 , 

n^oo |s|=n J-J- ' 3 
3 = 1 

where 7 is the speed of the additive branching random wedk with displacements In Y v , v G U^ =l {1, 2} n . 
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Proof. Simply write ]T'=i Y s\j = exp{£" In Yf,}. Then 


lim min (IT YAO 1 /” = exp{ lim min 

n—»oo \s\=n n—»oo |.s|=rc 

3 = 1 




n 


}• 


The assertion follows from the (BKH) theory since the exponential function is a continuous Injec¬ 
tion. □ 


This now provides the following approach to proving non-explosion by exploiting the theory 
of the speed of extremal (leftmost) particles in branching random walks. Namely, 

Proposition 5.3. For s G U^jl, 2} n , let { T s } be i.i.d. mean one exponentially distributed ran¬ 
dom variables independent of random variables in M 3 , and independent of { W s }. Assume that 

n 

lim inf min ( TT |W s i ? j _1 ) > 0. (5.5) 

n—>oo |s|=n ' ~~ 1 ’ 


Then P([£ = oo]) = 0. 

Proof In view of the Borel-Cantelli lemma one has that 


^ Pefrnin ^ |W s |j| 2 T s \j < M) < oo, for each M > 0 


. |s.|=n . 

n=l J=1 

is a sufficient condition for explosion to be a null event. Observe that for arbitrary M, A > 0, 


(5.6) 


P^minV |W sb |- 2 T s u < M ) = P( e - Amin h i=n£”=i |w sU | 2 Pu > e -^) 

I clrrrT) • ^ 

< e AM Ee _Amin i s l= n ^-'5= 1 l w sbP 2 Pb'. (5.7) 


|s|=r). 4- 
3 =1 


Also, since the mean of squares is larger than the square of the mean, and since the arithmetic 
mean is larger than the geometric mean, one has 


min^|W sb -| 2 T s \j = min^(|W sb -f 1 v / ^) 

Icl—n — |s|=n 


|a|=n 


3 =1 


3 =1 


|s|=n Ti * — 

3 = 1 
n 

> n min ( JJ |'W^l" 1 y/T ^) 2/ ”. 

s=n -*-4 

3 = 1 


(5.8) 


So the problem is reduced to showing that 


lim inf min ( TT iW^j 1 s/Tf\j) ” > 0. 

n-> 00 U V 


(5.9) 


3 =1 


25 



The indicated (positive) lower bound is possibly infinite. Since 


n 


n 


n 


min 

|s|=n 


( n i w *u 


-VtSS)> 


min 

|s|=n 


( n i w *u 

3 =1 


x ) pf in ( n v^b) 

s=n 

3 ! 


the two multiplicative factors can be treated separately. Moreover, the factor of n may in (5.8) 
may be included in either of these factors. The next calculation shows that it is most effectively 
assigned to the first factor. 

Namely, let 71 be the speed for JI"=i y/T s \j. Then 71 is directly computable from the above 
variant Proposition 5.2 on (BKH). However it is sufficient to bound 71 away from zero. Accord¬ 
ingly one has the following simple estimate. 

For M > 0 and u > 0 one has 


E II < M > ^ E^OIv^^) 

n= 1 j =1 


, I*l=re . . 

n= 1 J=1 


E 2 ”-P( Yn E ln ( T ‘w) < ln M ■ 


n— 1 
00 


3 = 1 


y 2 n P( e - ^ E)”=! ln ( T ib) > e ~ un in m ^ 


n= 1 
OO 

s E e 

n= 1 


n(ln 2 +m In M+ln T(l- ^)) 


(5.10) 


This series converges for T(1 — |) < \M u . Thus, taking u — 1, the series converges for any 
M < -V. It now follows from the Borel-Cantelli lemma that 

Ay/ 7T 


71 > 



> 0. 


n 

Regardless of the approach taken, the resolution of the explosion problem clearly involves 
a thorough understanding of the dilogarithmic branching random walk and its properties. We 
conclude this appendix with a few properties that may prove useful to this end and, at least, provide 
some insight into the technical nature of the problem in a future analysis. 

To the best of our knowledge, the dilogarithmic random walk is introduced in the present article 
for the first time. However an extensive treatment of the dilogarithm function, its properties and a 
selection of other applications in both physics and mathematics, is available in [19]. 

For the purposes of this article, let us note the invariance (multiplicative group symmetry about 
the identity) of the distribution of the ratios, one has 

p , l W *|j+il <r)= ( 27r- 2 [Li 2 (r) - Li 2 (— r)] if 0 < r < 1 

5 l w s|jl ~ \l - 27T- 2 [Li 2 (f) - Li 2 (-£)], if r > 1 . 

On the other hand, it is a rather direct calculation to check that 
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Proposition 5.4. 


E^ 
E^ In 
E^ (In 


|W Ail | 
|W S | 0 | 
|W A nl 
|W A ,|o| 
W sjl 


= oo, V 0 ^ 6 K 3 

= 0, VO/^GR 3 


W s | 0 | 


)' n < oo, Vm >1, 0 ^ G 


In particular, In |W s iJ = In |£| + V'” =l In n = 1 , 2 , , is a martingale. 

1 J 1 \ VV s |j_ 11 

As a consequence one has the following 

Corollary 5.1. The dilogarithmic random walk is 1-neighborhood recurrent in the sense that for 
fixed but arbitrary s G {1, 2 }°°,for each S > 1 

-P(|W A , n | <1 + 5 i.o.) = 1 

In particular, the path-wise explosion times are a.s. infinite for each path s. 

Proof By the Chung-Fuchs theorem it follows that In |W s i n | is O-neighborhood recurrent. That 
is, given e > 0 , 

P(| in |W A | n | < e i.o.) = P(e~ £ < |W S |„| < e e i.o.) = 1. 

The assertion follows by taking e = ln(l + 5). □ 

Corollary 5.2. 

a 2 R 2 2 a 

E - — -- = — arctan(^—), 6 > 0, a G M. 

a 2 R 2 + 9 n W6 

Proof. Define, g(x) = E r > ■ Justify differentiation under the integral to get, with c = 2/tt 2 , 


g\x) = cx 


—2 r 


in 


1 + r 


f 0 (r 2 + x 2 ) 2 1 — r 

= cx j w ((r 2 + a;2 r 1 ) in 


= cx lim 

£—>■0 + 

= cx lim 


fl-e 


Uo 


lim 

M—>• oo 


| dr 

1 + r 
1 — r 
M ^ d 


| dr 


l+e. 


Tk (( r 2 + X 2 )- 1 ) ln 

dr 


1 + r 


| dr 


-cx 


( lim 

ye-s> 0 + 


e->0+ \(1 — e ) 2 + X 2 
1-e i 


ln( 


2 — e N 


'0 


+ 


r 2 + a : 2 \ 1 + r 1 — r 


e (1 + e ) 2 + x 2 
1 

dr + lim lim 

6—y0+ Mwoo 


1 — r 
ln ( 2 + € ' 


e 

-M 


= —C 


lim 


1 + X 2 eW+ \ 

2 

lim lim 


(^arctan (r/x) + (x/ 2 )(ln( 


1 + r , 
1 — r' 


l+e 
r=1—e 


r 2 + x 2 + r r — 1 


dr 


r =0 


1 + X 2 e—►()+ A/ >oc y 

2 7T 2 1 


(^arctan (r/x) + (x/ 2 )(ln( 


1 + r , 

r — 1 ' 


M 

r=l+e 


= —C 


1 + x 2 2 n 1 T x 2 


(5.12) 
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Note that g{ 1) = 1/2. Indeed, one has 


1 _|_ r 2 


In 


r + 1 


r — 1 


dr 

r 


fOO r 2 


1 _|_ r 2 


In 


r + 1 


r — 1 


Thus, with c = 2 / 7 T 2 , 


3(1) = c 


5 1 + r 2 

r-oo r 2 


In 


r + 1 


= c 


In 


r — 1 
r + 1 


= cj In 


i 1 + r ; 

00 r + 1 


r — 1 
dr 

- = 1 / 2 . 

r 


dr f°° 1 , 

— + c / ——-2 In 

r J 1 1 + 

dr f°° 1 , 

— + c / —--2 In 

r ,/i 1 + r 


dr 

r 


r + 1 


r — 1 
r + 1 


r — 1 


dr 

r 

dr 

r 


r — 1 

Then, from (5.12) and (5.13) one has 

g(x) — — ( -arctanx') = — arctan(l/x) 

7r V 2 / 7r 

and the result follows setting x = Ve/a. 

Corollary 5.3. 

71 rr./ iri 

lim TT — 1 lt=1 —-— exists , 6 > 0. 

n ^f = \o+nuR? 

Proof. The limit exists by virtue of being a positive super-martingale. 


(5.13) 


□ 
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